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The recent WMAP data have confirmed that exotic daric matter together with the vacuum en- 
ergy (cosmological constant) dominate in the flat Universe. Thus the direct darlc matter detection, 
consisting of detecting the recoiling nucleus, is central to particle physics and cosmology. Supersym- 
metry provides a natural dark matter candidate, the lightest supersymmetric particle (LSP). The 
relevant cross sections arise out of two mechanisms: i) The coherent mode, due to the scalar inter- 
action and ii) The spin contribution arising from the axial current. In this paper we will focus on 
the spin contribution, which is expected to dominate for light targets. For both modes it is possible 
f~^ , to obtain detectable rates, but in most models the expected rates are much lower than the present 

f^ . experimental goals. So one should exploit two characteristic signatures of the reaction, namely the 

^\] ' modulation effect and, in directional experiments, the correlation of the event rates with the sun's 

motion. In standard non directional experiments the modulation is small, less than two per cent. In 
the case of the directional event rates we like to suggest that the experiments exploit two features, 
of the process, which are essentially independent of the SUSY model employed, namely: 1) The 
forward-backward asymmetry, with respect to the sun's direction of motion, is very large and 2) 
The modulation is much larger, especially if the observation is made in a plane perpendicular to the 
sun's velocity. In this case the difference between maximum and minimum can be larger than 40 
per cent and the phase of the Earth at the maximum is direction dependent. 
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PACS numbers: 95.35.+d, 12.60.J\ 



INTRODUCTION 



o 

^ ■ The combined MAXIMA-1 Q, BOOMERANG 0, DASI 0, COBE/DMR Cosmic Microwave Background (CMB) 

^^ , observations [^ , the recent WMAP data ,6] and SDSS |3 imply that the Universe is flat |5( and and that most of the 

^S^ '■ matter in the Universe is dark, i.e. exotic. 
Oh, 

X- nt = 0.044 ± 0.04, n,n = 0.27 ± 0.04, r^A = o.69 ± 0.08 

(D : . ..... 

r^ , for baryonic matter , cold dark matter and dark energy respectively. An analysis of a combination of SDSS and 
WMAP data yields [7| fim « 0.30 ± 0.04(1(7). Crudely speaking and easy to remember 

i>^ : 1^6 « 0.05, 17cDM» 0.30, 17a » 0.65 

5t 1 Since the non exotic component cannot exceed 40% of the CDM j^], there is room for exotic WIMP's (Weakly 

Interacting Massive Particles). In fact the DAMA experiment [9j has claimed the observation of one signal in direct 
detection of a WIMP, which with better statistics has subsequently been interpreted as a modulation signal [lOj . 

Supersymmetry naturally provides candidates for the dark matter constituents [ll|,[l3-ll3- ^^ ^^^ most favored 
scenario of supersymmetry the LSP can be simply described as a Majorana fermion, a linear combination of the 
neutral components of the gauginos and higgsinos lll|,ll3-llS|- ^^ most calculations the neutralino is assumed to be 
primarily a gaugino, usually a bino. Models which predict a substantial fraction of higgsino lead to a relatively large 
spin induced cross section due to the Z-exchange. Such models have been less popular, since they tend to violate the 
relic abundance constraint. The upper bound on this constraint has, however, been decreased by the recent WMAP 
data. In fact the LSP relic abundance (including co-annihilation) is: 

• Before WMAP: 

0.09 < ^Lsph^ < 0.22 

• After WMAP: 

0.09 < QLSph^ < 0.124 



These fairly stringent constrains, however, apply only in the thermal production mechanism. Furthermore they 
do not affect the LSP density in our vic inity derived from the rotational curves. We thus feel free to explore the 
consequences of two recent models |21|, [23, which arc non-universal gaugino mass models and give rise to large 
higgsino components. 

THE ESSENTIAL THEORETICAL INGREDIENTS OF DIRECT DETECTION. 

Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the observed rotational 
curves, it is essential to directly detect |ll|,[l3|-l2a| such matter. Such a direct detection, among other things, may 
also unravel the nature of the constituents of dark matter. The possibility of such detection, however, depends on the 
nature of its constituents. Here we will assume that such a constituent is the lightest supersymmetric particle or LSP. 
Since this particle is expected to be very massive, my. > 30GeV, and extremely non relativistic with average kinetic 
energy T « 50KeV{m^/l00GeV), it can be directly detected [ll|-[26j mainly via the recoiling of a nucleus (A,Z) in 
elastic scattering. The event rate for such a process can be computed from the following ingredients: 

1. An effective Lagrangian at the elementary particle (quark) level obtained in the framework of supersymmetry 
as described , e.g., in Refs |19ll2i^ . 

2. A well defined procedure for transforming the amplitude obtained using the previous effective Lagrangian from 
the quark to the nucleon level, i.e. a quark model for the nucleon. This step is not trivial, since the obtained 
results depend crucially on the content of the nucleon in quarks other than u and d. This is particularly true for 
the scalar couplings, which are proportional to the quark masses |27J — |8f]l| as well as the isoscalar axial coupling. 

3. Nuclear matrix elements |3l| — [33J obtained with as reliable as possible many body nuclear wave functions. 
Fortunately in the most studied case of the scalar coupling the situation is quite simple, since then one needs 
only the nuclear form factor. Some progress has also been made in obtaining reliable static spin matrix elements 
and spin response functions [23 ■ 

Since the obtained rates are very low, one would like to be able to exploit the modulation of the event rates due to the 
earth's revolution around the sun 34, 35] 361— 39]. In order to accomplish this one adopts a folding procedure, i.e one 
has to assume some velocity distribution ^, .Mji ^iii£|i[^-|^ for the LSP. In addition one would like to exploit 
the signatures expected to show up in directional experiments, by observing the nucleus in a certain direction. Since 
the sun is moving with relatively high velocity with respect to the center of the galaxy, one expects strong correlation 
of such observations with the motion of the sun |16l |44| . On top of this one expects to see a more interesting pattern 
of modulation as well. 

The calculation of this cross section has become pretty standard. One starts with representative input in the 
restricted SUSY parameter space as described in the literature for the scalar interaction IT, 19] (see also Arnowitt 
and Dutta |24j ) . We will only outline here some features entering the spin contribution. The spin contribution comes 
mainly via the Z-exchange diagram, in which case the amplitude is proportional to Z| — Z| (Z3, Z4 are the Higgsino 
components in the neutralino). Thus in order to get a substantial contribution the two higgsino components should 
be large and different from each other. Normally the allowed parameter space is constrained so that the neutralino 
(LSP) is primarily gaugino, to allow neutralino relic abundance in the allowed WMAP region mentioned above. Thus 
one cannot take advantage of the small Z mass to obtain large rates. Models with higgsino-like LSP are possible, but 
then, as we have mentioned, the LSP annihilation cross section gets enhanced and the relic abundance fl^ h? gets 
below the allowed limit. It has recently been shown, however, that in the hyperbolic branch of the allowed parameter 
space j2i|, |2l| even with a higgsino like neutralino the WMAP relic abundance constraint can be respected. So, even 
though the issue may not be satisfactorily settled, we feel that it is worth exploiting the spin cross section in the 
direct neutralino detection, since, among other things, it may populate excited nuclear states, if they happen to be 
so low in energy that they become accessible to the low energy neutralinos. In order to get simple estimates of the 
spin induced neutralino- nucleon cross-section under favorable circumstances, to be used as a guide to other slower 
processes (the directional rates, which is the main purpose of the present paper, and transition rates to the excited 
states to be studied elsewhere) we will utilize two recently proposed models: 

1. Non-Universal Gaugino mass models |21|. The gauginos belong to the adjoint of SU{5) , while the Chiral 
superfields $ are in the n-dimensional 
(24 X 24)sy™ = S n, n = 1, 24, 75, 200. 
The three gaugino masses at GUT are expressed in terms a single SUSY breaking parameter AI1/2, i-e. M- '" = 
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2 

M2 = -// 
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M2 = 2^^ 


3 


, Af2 = 3 



Cf Ml/2 The coefBcients Cf at the GUT scale are given in the table Q] One, of course, must take into account 
the renormalization effects. These are also given in the same table. The parameter fj, is constrained by: 
/i2 + 0.5Af| = r(")Af2^2 with r(") = (2.1, 0.3, 1.4) ^ n = (1, 75, 100) 

2. Anomaly mediated SUSY breaking (AMSB) |22|, inspired, e.g by superstring models [43 ■ Three simple cases 
will be considered: 

ABSM {Bino) : Mi = 
ABSM (Wino) : Mi = 

ABSM (Higgsino) : Mi 

with /i constrained by phenomenology alone. The linear relation between the Higgsino and the gaugino mass 
may not be a good approximation, furthermore it is not so easy to satisfy the LSP relic abundance constraint. 
We included, however, such a model in our calculation to compare its predictions on the neutralino direct 
detection rate with other models. 

In a given SUSY model one can calculate the amplitudes associated with the scalar and the spin amplitudes at the 
quark level. Going from the quark to the nucleon level, however, is not trivial. One distinguishes two cases: 

• The isoscalar axial current 

Most of the proton spin is not due to the quark spins (proton spin crisis-EMC effect). In fact one finds that the 
isoscalar and isovector axial current couplings transform as follows; 

/1('Z)-/1-3U^(9) ,/A(9)-/A-5UA(g) ,5^^0.1 , 5^ = 1-23 

Where the label (q) defines the quantity at the quark level. The axial current components /^, j\ are defined 
in the standard weak interaction formalism, see e.g. [23j. In other words they have been normalized so that 
CTp^^o — 3(/^ + F4)^cro, ffi^t^o = 3(/^ — i^j^)^cro, fo = j^^iGprrip)'^ , for the proton and neutron spin cross sections 
respectively. 

• The scalar amplitude 

The relevant amplitude at the quark level is proportional to the quark mass. Thus the naive quark model, in which 
one considers only u and d quarks in the nucleon, is not even approximately good. In going to the nucleon level one 
must compute the matrix element of the quark number operator multiplied by its mass, which we express as a fraction 
of the proton mass, i.e. : 

NlggmglN^^ = fgiTiN (1) 

Since we are primarily interested in the spin contribut ion, we will not elaborate here on how one obtains the quantities 
fq , but we wiU refer the reader to the literature 113, [Ij, Hi], |i3i El- 

Once the LSP-nucleon cross section is known, the LSP-nucleus cross section can be obtained. The differential cross 
section with respect to the energy transfer Q for a given LSP velocity v can be cast in the form 



du 
2{JIrbvY 



d(T{u,v) = -^-^-—[{tsF\u) + t,p,nFiiiu)] (2) 



where we have used a dimensionless variable u, proportional to Q, which is found convenient for handling the nuclear 
form factor [26| F(u), namely 

u^-^ , Qo « 40 X A-^^^ MeV. (3) 

Qo 

Hr is the reduced LSP-nucleus mass and b is (the harmonic oscillator) nuclear size parameter. In the above expression 
we have neglected the small vector and pseudoscalar terms. 



Furthermore 

V.c = rr^ „ 4^ 

firip) ~ rup is the LSP-nucleon reduced mass and 



//2 



M^ \2 spin 
>, TT/ ^« v*^ -^ spin -) ^ spin — ^o 

^ir{p) P'^ '^ '^ 3(1 ,4^2 



^spin — \ f \) ^T) v*^ sspinTSspin — ^o 0\U) \p) 

" "" ' 3(1 + ^)" 

J A 



CT and (T^ are the proton cross-sections associated with the spin and the scalar interactions respectively and 

Sin) = [(:^f7o(0))^|4l + 2^"o(0)17i(0)|4^ + nm? ] (6) 

The precise definition of the spin response functions Fij, with i,j — 0,1 isospin indices, which are essentially the 
"spin form factors" normalized to unity at zero momentum transfer, can be found elsewhere [33J. As we have already 
mentioned the existing experimental limits imply that the scalar LSP-nuclcon cross section satisfies: af^ q < 10~^pb. 
The constraint on the corresponding spin cross-section is less stringent. 

Some static spin matrix elements J23 , ^M ' |2a| for some nuclei of interest are given in table Hll 

RATES 

The differential (non directional) rate with respect to the energy transfer u can be written as: 

dRundtr^ z da{u,v)\v\ (7) 

Where p{0) = 0.3GeV/cm^ is the LSP density in our vicinity, m is the detector mass, rriy^ is the LSP mass and da{u, v) 
was given above. 

The corresponding directional differential rate, i.e. when only recoiling nuclei with non zero velocity in the direction 
e are observed, is given by : 

dRdtr = -. \v\v.e Q{v.e) -— da{u,v) (8) 

m^ ArriN 2,1: 

The LSP is characterized by a velocity distribution. For a given velocity distribution i{v'), with respect to the 
center of the galaxy. One can find the velocity distribution in the lab frame f{v, Ve) by writing 
v'=v+ Ve , Ve=Vq+Vi 

Vq is the sun's velocity (around the center of the galaxy), which coincides with the parameter of the Maxwellian 
distribution, and Vi the Earth's velocity (around the sun). The velocity of the earth is given by 

Ve = Vqz + 1^1 ( sina x — cosa cos^ y -I- cosa siwj z ) (9) 

In the above formula z is in the direction of the sun's motion, x is in the radial direction out of the galaxy, y is 
perpendicular in the plane of the galaxy {y = z x x) and 7 « 7r/6 is the inclination of the axis of the ecliptic with 
respect to the plane of the galaxy, a is the phase of the Earth in its motion around the sun (a = around June 2nd). 
The above expressions for the rates must be folded with the LSP velocity distribution. Wc will distinguish two 
possibilities: 

1. The direction of the recoiling nucleus is not observed. 
The non-directional differential rate is now given by: 

/ dRundir \^ _ /^\ _ P(0) "t /T^('^ 

where 



du I \ du I mv Am^ du 



,dT, f \v\ da{u,v) , 

' \- ' --f{'v,VE) \ ' d^v (11) 



du J \/JiP) ' du 



2. The direction e of the recoihng nucleus is observed. 
In this case the directional differential rate is given by: 

,,dE, > f v.e Q(v.e) „, .da(u,v) ,^„, 

{{-r)d^r) = j=^f{v,VE) ^ ' (12) 

-— d(-^^- v.ejd'^v 

2tt firbv 

The above coordinate system, properly taking into account the motion of the sun and the geometry of the galaxy, is 
not the most convenient for performing the needed integrations in the case of the directional expressions. For this 
purpose we go to another coordinate system in which the polar axis, Z, is in the direction of observation (direction 
of the recoiling nucleus) via the transformation: 

X \ / cos 6 cos $ cos 8 sin $ — sin 8 

F == -sm<i> cos$ 

Z J \ sin 8 cos $ sin 8 sin <i> cos 8 
It is thus straightforward to go to polar coordinates {6 , (p) of the new system in velocity space and get: 




\du/dir m^ AniM JavoV^ Jo Jo 27r firbv du 

with ^ — cos9 and A is the asymmetry parameter, which enters when one makes the replacement u^ =^ (1 + A)w|, i = 
y, z in the Maxwell-Boltzmann distribution (with respect to the galactic frame), while the x-component remains 
unchanged, see e.g. (38J. The orientation parameters 8 and <i> appear now explicitly in the distribution function / 
and not implicitly via the limits of integration. The function / can be obtained from the velocity distribution in a 
straight forward fashion. The (5 function ensures that in the directional case the variables u, v and ^ obey the required 
relation. In our numerical calculation we found it more convenient to use y — v Jvq and to express ^ in terms of 
the other two, namely u and y. So one is left with two integrations, over (/> and y. This way, to leading order in 
S = 2^ = 0.27, we find: 

/(A, 8, $, e, 0, V, ve)S{^ - ^ ff(A, 8, $, 0, aV^, y, S) (14) 

with 

g{\,Q,^,(j),z,y,S) = E{e,^,X,(f>,z,y)[l + S sina ^(8, $, 0, z, y)) + 5 cos a B{Q,^, X,(f>,z,y)] (15) 

where z — a^/u and a is given by: 

a = [firbvoV2]^ 
Furthermore 

S(8,$,A,0,z,y) - (l + A)[-sin7- 

sin7(zcos8 — \/y'^ — z^ cos <j) sin 8) + cos 7 (z sin 8 sin $ + 

^y2 „ z^ (cos $ sin + cos8cos0sin$))] (16) 



j4(8, ^,(l),z,y) — — zcos$sin8 — y j/^ — z^ (cos 8 cos cos $ — sin0sin$) (17) 



-(l + A) l+a^+2 z cos e -yV-z^ cos sine 

EiQ,<i>,(j),z,y)^e V V V/ x 



A I —z cos ^ sin Q — -w y'^ — z^ {cos cos (p cos $— sin (p sin $) 

e V / (18) 

To obtain the total rates one must integrate the expressions II 1 1 and 1 1 31 over the energy transfer from Qmin determined 
by the detector energy cutoff to Qmax determined by the maximum LSP velocity (escape velocity, put in by hand in 
the Maxwellian distribution), i.e. 

Vesc = 2.84 vo , vo = 229 Km/s. (19) 

In our analysis we included only the rotational velocity of the sun around the center of the galaxy. The local component 
of the sun's velocity is 10 times smaller and it can be neglected. 
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FIG. 1: Scatter plots associated with the scalar interaction with non-universal parameters obtained Cerdeno, Gabrielli, Gomez 
and Munoz. The light grey dotted area corresponds to the laboratory constraints, while the dark dotted area is associated with 
old relic abundance constraint and the black dotted area with the WMAP relic abundance constraint. 
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FIG. 2: The Total detection rate per {kg — target)yr vs the LSP mass in GeV for a typical solution in our parameter space 
in the case of ^^^/. Thick lines correspond to model B, while fine lines to model C. In the upper curve no detector cutoff was 
employed, while in the lower curve we used a detector energy cutoff of Qmin = 10 KeV . Such effects introduce variations in 
the rates by factors of about two. 



RESULTS 



We will specialize the above results in the following cases: 
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FIG. 3: The quantity R, ~ event rate for Qmin = 0, associated with the spin contribution in the case of the A — 19 system 
for model n=l of Chattopadhyay and Roy on the left and model n=75 on the right. The curves from top to bottom correspond 
to {sign{jj,) , tan{f3)) : (— , 40), (+, 40), (— , 10), (+, 10). The curves for ±10 cannot be distinguished from each other. 
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FIG. 4: On the left the same as in Fig. Inland on the right the proton cross section in pb, both for n — 100. Since the isoscalar 
contribution is negligible, the neutron cross section is similar. 



Non directional unmodulated rates 



Ignoring the motion of the Earth the total non directional rate is given by 



R = Rt{a,Qmin) , R — 



P(0) 



m^Q Anip 



flr{p) 



f VP>[ 



A' 



spin 
To 



sptn \ 



(20) 



where t is the ratio of the calculated rate divided by that obtained using the previous equation. It represents the 
modification of the total rate due to the folding procedure and the nuclear structure effects, t depends on Qmim i-e. 
the energy transfer cutoff imposed by the detector and the parameter a introduced above. All SUSY parameters, 
except the LSP mass, have been absorbed in R. 

Via Eq. H2U|) we can, if we wish, extract the nucleon cross section from the data. For most of the allowed parameter 
space the obtained results for the coherent mode are undetectable in the current experiments. As it has already been 
mentioned it is possible to obtain detectable rates. Such results have, e.g. recently been obtained by Cerdeno et al 
[49j with non universal set of parameters and the Florida group [50| . 

A representative set is shown in Fig. Qor with universal couplings 17] for large tan/3 in Fig. I21in the case of the 
targeti27/. The planned experiments, Hke CDMS 51], EDELWEISS _52J, IGEX 53], ZEPLIN _54] and GENIUS 54 
are however, expected to improve so that they may detect rates two or three orders of magnitude smaller (see Fig. ^). 

In the case of the spin contribution we see that detectable rates are possible only in those cases in which the 
neutralino has a substantial Higgsino component (see Figs OJ- 0] and Figs EJ- Ell . Our results for the nucleon cross 
section, obtained by neglecting the isoscalar axial current due to the EMC effect, are in essential agreement with more 
detailed calculations, which have appeared after our manuscript was completed 20]. 

For the target ^^^/ the corresponding results are shown in Figs [71- |S1 and Figs El- E| in the case of the spherically 
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FIG. 5: The quantity R, ~ event rate for Qmin = 0, associated with the spin contribution in the case of the A = 19 system 
for model Mi = (l/3)/i, M2 = (2/3)^ of Murakami and Wells on the left and model Mi = (3/2)^, M2 = (1/2)^ on the right. 
The curves are independent of the sign of /i. From top to bottom correspond to tan{f3) : 40, 10. 
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FIG. 6: The same as in Fig. Ofor model Mi — (3/2)/i, M2 — 3/i of Murakami and Wells. Note the difference in the scale of the 
LSP mass compared to Fig. |3 



symmetric M.B. velocity distribution and Qmin — 0. The effects of asymmetry and the detector energy cutoff Qr, 
are analogous to those previously found in the coherent mechanism. 
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FIG. 7: The event rate, associated with the spin contribution in the case of the A = 127 system. The notation is the same as 
in Fig. 121 Only the spherically symmetric velocity distribution (A — 0) has been considered. 



Rspin 



4- 



2 




100 120 140 160 180 200 



mlsp-> 



FIG. 8: The same as in Fig. |7|for model n=100. 
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FIG. 9: The event rate, associated with the spin contribution in the case of the A — 127 system for modelMi = (l/3)/i, M2 = 
{2/3)fi of Murakami and Wells on the left and model Afi = (3/2)/i, M2 = (1/2)(U on the right. The curves are independent of 
the sign of /i. From top to bottom correspond to ian(/3) : 40, 10. 



Modulated Rates 

If the effects of the motion of the Earth around the sun are included, the total non directional rate is given by 

R = Rt[{l + hia,Qrmn)cosa)] (21) 

with h the modulation amplitude, relative to the unmodulated (time averaged) amplitude, and a is the phase of the 



Rspin 




100 120 140 16 



180 200 



mlsp-> 



FIG. 10: The same as in Fig. |S]for model Mi = (3/2)^, M2 = 3/i of Murakami and Wells. Note the difference in the scale of 
the LSP mass compared to Fig. |5| 
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FIG. 11: The parameter tdir defined in the text as a function of the angle B for A = 19 at the top and A — 127 at the bottom. 
On the left we show tdir for A = and on the right for A = 1. The parameter tdir is independent of <3? for A = and it depends 
slightly on it, when A is different from zero. The results presented correspond to an LSP mass of 100 GeV. 



Earth, which is zero around June 2nd. The modulation amplitude would be an excellent signal in discriminating 
against background, but unfortunately it is very small, less than two per cent (see table Ullll . Furthermore for 
intermediate and heavy nuclei, it can even change sign for sufficiently heavy LSP. So in our opinion a better signature 
is provided by directional experiments, which measure the direction of the recoiling nucleus. 



Directional Rates. 



Since the sun is moving around the galaxy in a directional experiment, i.e. one in which the direction of the 
recoiling nucleus is observed, one expects a strong correlation of the event rate with the motion of the sun. In fact 
the directional rate can be written as: 



trl - K - 

Rdir = -7r~R [1 + h,nCos{a - am tt)] = --R t[l + h„ 
Zn Ztt 



cosya — a,-, 



(22) 



where tdir is a quantity analogous to t discussed above, and /i,„ is the modulation, am is the "shift" in the phase 
of the Earth a, since now we have both sine and cosine terms and the maximum occurs aX a — amTr. k/{2t:) is the 
reduction factor of the unmodulated directional rate relative to the non-directional one. The parameters n , hm , Q^m 
depend on the direction of observation: 

e = (sin O cos $ , sin 8 sin <& , cos 9) 

The parameter tdir for a typical LSP mass 100 GeV is shown in Fig. E]as a function of the angle 8 for the targets 
A = 19 and A = 127. We see that the change of the rate as a function of the angle 8 for the Maxwellian LSP velocity 
distribution is quite dramatic. This figure is important in the analysis of the angular correlations, since, among other 
things, there is always un uncertainty in the determination of the angle in a directional experiment. 

We sometimes prefer to use the parameters k = ^^ and hm, since, being ratios, are expected to be less dependent on 
the parameters of the theory. We first exhibit the dependence of hm on the angle 8 for an LSP mass of m^ = lOOGeV 
in Figs 1121 and 1131 Then we exhibit the dependence of the parameters i, /i, K,hm, and a™, which are essentially 
independent of the LSP mass for target A — 19, in Table IIIII (for the other light systems the results are almost 
identical). 

The asymmetry is quite large. For a Gaussian velocity distribution we find: 

A.= ^i-i-^;+-i .o.97 

R{-z) + R{+z) 
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FIG. 12: The modulation amplitude h,n in the case of the A = 19 system as function of O for <1> = 0, 27r at the top, $ — -k 12 
in the middle and <l> = 3 7r/2 at the bottom. The values for angles < 0.5 should be discarded since the rate is tiny (see Fig. 
IIH . Otherwise the notation is the same as in Fig llll 



In the other directions it depends on the phase of the Earth and is equal to almost twice the modulation. For a 
heavier nucleus the situation is a bit complicated. Now the parameters k and hm depend on the LSP mass. The 
situation is exhibited in Figsll4landll5l The asymmetry and the shift in the phase of the Earth are similar to those 
of the A—V^ system. 

In the results shown we only considered the rotational velocity of the sun around the center of the galaxy. We do 
not, however, expect a significant modification, if the other components of the sun's velocity, which are an order of 
magnitude smaller, are included. 



CONCLUSIONS 



It is well known that, in the coherent case, only in a small segment of the allowed parameter space the rates 
are above the present experimental goals llTl 113 . |2J| , which of course may be improved by two or three orders of 
magnitude in the planned experiments |5l|-|55|. In the case of the spin contribution only in models with large higgsino 
components of the LSP one can obtain rates, which may be presently detectable, but in this case, except in special 
models, the bound on the relic LSP abundance may be difficult to respect. Anyway it appears that in both cases the 
expected rates are small. It may, therefore, be necessary to exploit any characteristic experimental signatures, which 
may reduce the formidable backgrounds at such low counting rates. 

In the present paper we considered the spin induced rates and examined the following signatures: 

• Correlation of the event rates with the motion of the Earth (modulation effect) 

• Angular correlation of the directional rates with the direction of motion of the sun as well as their seasonal 
variation. 



12 



T 

g 






e 



3 e- 



T 



1" 0.15 

J 0.1 ■ 




e 




0.5 1 



3 e. 




e 



FIG. 13: The same as in Fig. [Hlfor A=127. 
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FIG. 14: The parameter k as a function of the LSP mass in the case of the A — 127 system, associated with the Maxwellian 
velocity distribution (A — 0) and Qmin ~ 0. On the left we show the results for O — -k 12 and on the right = tt (the results 
are independent of $). The rate for O = is negligible. 



Such experiments are currently under way, like the UKDMC DRIFT PROJECT experiment |44 
Detector of the Kyoto- Tokyo collaboration 57] and the TOKYO experiment |58l |. 



the Micro-TPC 



Let us first consider the conventional experiments and focus on the relative parameters t = RjR and the modulation 
(seasonal variation) amplitude /i, normalized to zero when the motion of the Earth is ignored. In the case of light 
nuclear targets they are essentially independent of the LSP mass, but they depend on the energy cutoff, Qmin- For 
Qmin = they are exhibited in Table ITTTl They are essentially the same for both the coherent and the spin modes. 
For intermediate and heavy nuclei they depend on the LSP mass [56j- It is clear that among the light targets, from 
the point of view of the static spin matrix elements, the most favored system J33j is the A — 19, since the spin matrix 
elements are both large and very reliably calculated. We should keep in mind, however, that for heavy LSP the 
reduced mass can be large in the case of a heavy nucleus like ^°^P6. The increase of the rates caused by the increase 
of the reduced mass may very well compensate for the smallness of the spin matrix elements. As a matter of fact in 
the case of ^^^/ we see that, for an LSP mass greater than 100 GeV, the obtained rate is larger than that of ^^F in 
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FIG. 15: The modulation amplitude hm- From left to right and top to bottom the graphs correspond to ±x,+y,—y,—z. 
Otherwise the notation is the same as in Fig 1141 



spite of the fact that the dominant isovector static spin matrix elements are in the opposite order. We should not 
forget, though, that the calculated spin matrix elements in the case of the A — 19 system are more reliable. 

In the case of the directional experiments, i.e when only nuclei recoiling in a certain direction are counted, we can 
summarize our results as follows: 

• The angular dependence of the rate. 

The ratio of the directional rate divided by that of the usual rate , given by « k/{2t:), is smaller than I/Stt. 
Such a big loss in the rate may very well be compensated by the important experimental signature associated 
with the fact that the factors k depend on the direction of observation. In the case of a Gaussian velocity 
distribution this factor is the largest, n « 0.4, when the nucleus is recoiling opposite to the direction of motion 
of the sun. k is the smallest in the direction of the sun's motion (least favored direction) and in between for 
the other directions. In other words the directional rate is strongly peaked in the direction opposite to the 
velocity of the sun and the resulting event asymmetry between these two directions is very large. In a plane 
perpendicular to the sun's direction of motion k takes intermediate values, k « 0.08. The asymmetry between 
one direction and its opposite is in this case almost zero, if the motion of the Earth is neglected. Anyway the 
study of such angular correlations may play a role in confirming any observed neutralino events, since there may 
appear seasonal effects, which can mimic the small modulation in the conventional experiments. 

• Seasonal variation (modulation) of the directional rates. 

We have shown that the directional rates can, in addition, exhibit seasonal variation (modulation). In the most 
favored direction the modulation is not very large, but still it is three times larger compared to that expected 
in the standard non directional experiments. This gain is , perhaps, not big enough to compensate for the 
reduction in the number of the events. In a plane perpendicular to the sun's direction of motion, however, the 
modulation is quite large (see Table llll() . The experimentalists will themselves decide what use, if any, to make 
of these predictions. A simple argument indicates that they may be useful. Taking, as an example, events in 
the x-direction (radial galactic direction), we can see that the modulation signal will be reduced by a factor of 
(27r/0.08) (0.272/0.02) = 5.8 The background will be reduced by a factor of v/(27r/0.08) = 8.9. Thus one expects 
a gain of 1.5. 

Furthermore the modulation in this plane is characterized by a very interesting seasonal pattern, depending on 
the angle of observation (see Table Ulljl . Perhaps this also can be exploited by the experimentalists. 

The predicted reductions in the rates of the directional experiments compared to the standard experiments, at a 
first sight, may be seen as an obstacle to be overcome for the benefit of the above good signatures. This may be 
true in the case of some of the planned experiments |58l | , which intend to use organic detectors capable of making 
observations in only one predetermined direction for each run. Quite clearly, however, the directional observations 
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are going to offer only advantages in the case of experiments using TPC detectors |4J,|57|. The TPC counters can 
simultaneously register all events. If something interesting is found, the analysis can be made directionally to reject 
possible background events. 

We finally hope that, in spite of the reduction in the predicted rates, the signatures of the directional experiments 
(large asymmetry and modulation amplitude) can be exploited by the experimentalists. 
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TABLE I: A summary of the parameters C" of the model of Chattopadhyay and Roy and the relevant renormalization factors. 



n/i 


3 2 1 


comment 


1 

24 
75 
100 


1 1 1 
1 -3/2 -1/2 
1 3 -5 
1 2 10 


SUGRA (Bino) 

Bino 

Zino 

Higgsino 




25/9 25/30 25/60 


Ren. factor 



TABLE II: The static spin matrix elements for various nuclei. For light nuclei the calculations are from Divari et al (see text) 
. For ^^"^ I the results are from Ressel and Dean (see text) (*) and the Jyvaskyla-Ioannina collaboration (private communica- 
tion) (**). For ^°^ Pb they were obtained previously (see text). 



^Si ^^'Na ^^'r '^'r* 



'Pb+ 



[no(0)]2 2.610 


0.207 0.477 3.293 


1.488 


0.305 


[ni(0)]2 2.807 


0.219 0.346 1.220 
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0.231 


^(,(0)^1(0) 2.707 


-0.213 0.406 2.008 


1.501 


-0.266 


Mth 2.91 


-0.50 2.22 






Me.p 2.62 


-0.56 2.22 






A»ea:p 


0.99 0.57 
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TABLE III: The parameters t, h, k, hm and Qm for the isotropic Gaussian vefocity distribution and Qmin = 0. The results 
presented are associated with the spin contribution, but those for the coherent mode are similar. The results shown are for 
the light systems. For intermediate and heavy nuclei there is a dependence on the LSP mass. +x is radially out of the galaxy 
(0 = 7r/2, <1? — 0), +2: is in the sun's direction of motion (O = 0) and +3/ is vertical to the plane of the galaxy (O — 7r/2, $ — 7r/2) 
so that {x, y, z) is right-handed, am = 0, 1/2, 1, 3/2 means that the maximum occurs on the 2nd of June, September, December 
and March respectively. 

type t h dir k hm ctm 

+z 0.0068 0.227 1 

dir +(-)x 0.080 0.272 3/2(1) 

+(-)y 0.080 0.210 (1) 

-z 0.395 0.060 

aU 1.00 

aU 0.02 



